Analysis and interpretation of high transverse entanglement in optical parametric 

down conversion 
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Quantum entanglement associated with transverse wave vectors of down conversion photons is 
investigated based on the Schmidt decomposition method. We show that transverse entanglement 
involves two variables: orbital angular momentum and transverse frequency. We show that in 
the monochromatic limit high values of entanglement are closely controlled by a single parameter 
resulting from the competition between (transverse) momentum conservation and longitudinal phase 
matching. We examine the features of the Schmidt eigenmodes, and indicate how entanglement can 
be enhanced by suitable mode selection methods. 

PACS numbers: 



Nonclassical correlation among photons has been one 
of the central topics in quantum optics and quantum 
communication. In particular, the property of quantum 
entanglement is now recognized as the key to accomplish- 
ments exceeding the limit imposed by classical physical 
laws. Entangled photons are experimentally available 
through optical spontaneous parametric downconversion 
(PDC). Such highly correlated photons have been em- 
ployed in fundamental experiments, e.g., to test the viola- 
tion of Bell's inequality 1], and to demonstrate quantum 
teleportation P|. Significantly, down conversion pho- 
tons provide quantum entanglement in a rich diversity of 
forms beyond polarization Bell states. For example, one 
has continuous entanglement in field frequencies [j, 0, 
energy-time variables jsj and momenta 6]. Altogether, 
counting the discrete polarization and orbital angular 
momentum variables j^, J] , a full quantum characteri- 
zation of just two-photon states already requires multiple 
parameters and the term 'hyper- entanglement' has been 
proposed in order to emphasize the complex structure of 
such quantum states |inl |. 

Here we focus on PDC entanglement that is associ- 
ated with the transverse wavevectors of signal and idler 
photon pairs. We are especially interested in ways to 
reach states of anomalously high entanglement. The 
PDC quantum state is sometimes called a biphoton ■ 
Biphoton correlation is experimentally accessible and 
several investigations have addressed this topic in re- 
gard to the underlying orbital angular momentum struc- 
ture 0, 0, 01 and applications in quantum imaging [T^ 
are being explored, including comparison with two-beam 
correlation using spatially coherent classical light |Q. 
In this note we approach the subject using the single- 
sum Schmidt decomposition technique which yields 
for the first time a complete characterization of the ex- 
isting entanglement by determining the natural set of 
bi-orthogonal mode pairs. This is in contrast to stud- 
ies based on familiar double-sum expansions in Gauss- 



Laguerre or Gauss-Hermite modes. 

To begin with, we consider the signal and idler pho- 
tons of PDC under the so-called type-II phase matching 
condition. The two photons are distinguished by their 
orthogonal polarizations. The explicit forms of bipho- 
ton states in three dimensional space are complicated 
by the details of crystal properties and pump pulse pro- 
files 01 . However, under plausible assumptions in the 
paraxial approximation, the dependence of longitudinal 
components of the wavevectors kz and niay be elimi- 
nated in the monochromatic limit U ^3 ■ Denoting k_L 
and as transverse wavevectors of the signal and idler 
photons, the biphoton state takes the approximate form: 
I*) = / dkj_ / dq± C (k| ,q I ) |k | ,q I ) . Here the bipho- 
ton amplitude is given by |3, 113 , 



C(ki,qx)=£p(ki+qi)A(kj 



(1) 



where £p(k± + qj_) describes the transverse profile of the 
pump field in wavevector space, and A(kj^ — q^) is a 
purely geometrical function that will be specified later. 
Note that if the amplitude C(kj_,qj_) were to separate 
into factors depending only on k^ and q^ , then the state 
I^P) would be factorable (not entangled). 

The Schmidt decomposition of C (k^ , q^ ) corresponds 
to the expansion. 



A„u„(kj_)w„(qj 



(2) 



where Mn(kj_) and w„(q_L) are Schmidt modes defined 
by eigenvectors of the reduced density matrices for the 
signal and idler photons respectively, and A„ are the cor- 
responding eigenvalues |0| . Thus the mode functions u„ 
form a complete and orthonormal set, and the same is 
true for the w„ . Because density matrices always have fi- 
nite trace, the Schmidt decomposition is always discrete, 
even when the original specification of the state is natu- 
rally continuous, or doubly continuous as in the present 
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situation. This discreteness is intrinsic, independent of 
any artificial box boundary conditions. 

The Schmidt decomposition ^ provides the informa- 
tion eigenstates of the two-particle system, which char- 
acterizes the structure of entanglement by providing two 
important pieces of information not available otherwise. 
First, the Schmidt single-sum form of decomposition re- 
veals exactly how the photons are paired, i.e., if a signal 
photon is detected in the mode w„, then with certainty 
the idler photon must be in the mode «„. Second, the 
Schmidt eigenvalues A„ serve to measure the degree of 
entanglement. This is usually discussed in terms of the 
entanglement entropy E = — ^'-'£2 ^n- However, a 

more transparent and experimentally more direct mea- 
sure of entanglement is the 'average' number of Schmidt 
modes involved. The Schmidt number (or participation 
ratio) K provides this average: K = l/X^n'^n- The 
larger the value of K, the higher the entanglement [lij . 
We remark that the maximum value of K is governed by 
the volume of accessible phase space of the system under 
relevant physical constraints such as energy-momentum 
conservation. 

We first examine a class of biphoton amplitudes in 
which both £p and A are represented by gaussian func- 
tions: 

Cg(ki,qi) = AAge-l'^^+'i^l'/'^'e-^'I'^^^i^l'. (3) 

Here J\fg is a normalization constant, and the parameters 
a± and are the widths of the respective gaussians, to 
which we give physical meaning below. With this ampli- 
tude, it can be shown that the Schmidt eigenmodes are 
simply the energy eigenfunctions of a two-dimensional 
isotropic harmonic oscillator. The probability of finding 
a mode (in polar coordinates) with radial and angular 
quantum numbers n and m is proportional to ^2n+|T?x|^ 
where £, = {I — ba±)'^/ (1 -I- ba±)^. This leads to a very 
convenient closed form for the Schmidt number: 

Thus, K increases as the "control parameter" ba± in- 
creases or decreases, i.e., for both baj_ » 1 and ba± <C 1 
(see Fig. 1). It is interesting that K can be expressed 
in terms of the variance of the transver se wavevector 
components, K^^'^ = (g^) / \J where s±j = 

{kj ± qj)/\/2 with the subscript j referring to a com- 
ponent on the transverse xy plane. Such an expression 
connects the abstract measure of entanglement with ob- 
servable fluctuations. 

The biphoton amplitude Cg can be recognized as a con- 
tinuous representation of two-mode squeezed states, and 
these states maximize the EPR correlation for a given 
amount of entanglement Eq. (4) provides a use- 

ful estimation of degree of entanglement according to the 
widths of the functions £p(k^ + q^) and A(kx — q±)- 
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FIG. 1: Participation ratio X as a function of ba± for the 
biphoton amplitude Cg (dash line) and Ca (solid line). 

In practice, A is determined by phase matching in the 
longitudinal direction, and it is generally not a gaussian. 
We now follow Monken, et al. [ig and consider another 
class of biphoton amplitude: 

C,(k^,q^) -AA,e-l''-+i-l'/-lsinc [-b^\k^ - q^\'] . 

(5) 

Here the transverse profile of the pump field £p is treated 
as a gaussian, and this gives physical meaning to the 
width a±, and J\fs is a normalization constant. We 
have used the subscript s in order to emphasize that 
A(k^ — q^) is now taken as the sine function arising 
from longitudinal momentum mismatch in the PDC crys- 
tal. Spatial coherence is the physical effect that deter- 
mines b via the relation 6^ = cL/4ujp, where L is the 
length of the crystal and Up is the pump frequency. 

The amplitude © is an approximation to actual bipho- 
ton states realized in laboratories. This is because tech- 
nical details involving dispersion and birefringence effects 
are omitted. However, it captures the main features im- 
posed by conservation of energy and momentum, which 
are the key constraints on the state in the two-photon 
Hilbert space. To more fully appreciate © in terms of 
conservation rules, one can see first that the gaussian 
with the argument kj^ + qj^ is merely a statement of the 
uncertainty in transverse momentum conservation, and 
second that the sine function's argument expresses en- 
ergy conservation uig ^ Wi ~ LLjp/2 in practice p^ . 

In several previous studies of biphoton states, the sine 
function was approximated by unity in the small L limit 
(or equivalently the Fourier transform of A was approx- 
imated by a delta function) However, such an 
approximation is inapplicable to our Schmidt analysis 
here because the corresponding Cs{k±,q±) would be un- 
bounded on the kj^ = — q± manifold. Since the sine 
function ultimately limits the accessible transverse wave 
vectors, it has to be fully accounted for in the equation. 
We also note that exact analogs of the sum and difference 
arguments in jSJl appear in wave functions of entangled 
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states of two massive particles when expressed in terms 
of center of mass and relative coordinates ji^l ■ 

To carry out the complete Schmidt decom- 
position of (O, we make use of polar coordi- 
nates, writing: kj^ = (k±cos9k, fcj_sin0fc) and 
qj^ = {q± cos 9q, q±sm9q). The biphoton ampli- 
tude can be decomposed in the form: Cs(k^,qj_) = 
E™e™(^^-''')V^^F,„(A:^,9i)/2^, where F^{k^,q±) is 
a function of magnitude k± and q±, and is normalized 
according to \Fjn{k±, q±)\^k^qj^dk^dq±^ = 1. We 

identify the integer m as a quantum number of orbital 
angular momentum. Therefore P,n is the probability of 
finding the two photons with opposite orbital angular 
momentum numbers m and —m. 

Further decomposition of Fm{k±, q±)\/k±q± is a sim- 
pler task due to its lower dimensionality as com- 
pared with Cs(kj_,qj_). According to the Schmidt de- 
composition scheme, we have: F„i{k±, qi_)^/kj_q± = 
{k±)(j)n,m{q±), where the (f>n,m are the 
same for idler and signal photons because of the sym- 
metry assumed in Eq. jS)). Consequently the Schmidt 
decomposition of the biphoton amplitude reads: 

oc oc 

(qx) 

m— — oo n— 

(6) 

where Xnm = Pmlnm and Un,-m{^±) = 

e^"^^''(j)n^m{k±)/V^T^k± are normalized Schmidt mode 
functions in wavevector space. Eq. (6) reveals that 
quantum entanglement involves orbital angular momenta 
and the magnitude of transverse wavevectors. 

The gaussian and the sine function in (jSJ compete with 
each other by enforcing correlations in opposite direc- 
tions. High entanglement is achieved only when one of 
them becomes dominant, i.e., a±b ^ 1 or a±b ^ 1. We 



have numerically performed the Schmidt decomposition 
for various values of aj_b. In Fig. 1, we show the Schmidt 
number K ~ 1/ X^nm ^nm as a function a±b. Consistent 
with the description above, we see that either a decrease 
or increase in a±b can lead to an increased degree of en- 
tanglement. For example, K 7.2 when a±b = 0.3. The 
minimum Schmidt number K w 1.4 occurs at a±b w 1. 
The figure also indicates that Cs is more entangled than 
Cg at the same parameter value <T±b. 

To learn how various modes contribute to the entan- 
glement, we plot in Fig. 2 the distribution of for the 
case ba± — 0.25. Such a control parameter value corre- 
sponds to a moderately high entanglement regime where 
conservation of transverse momentum dominates. Small 
values of ba± are typically accessible in experiments us- 
ing very thin crystals. In Fig. 2 we clearly see that most 
probabilities are found on the n — manifold, where the 
quantum number n equals the number of nodes in the 
radial direction (Fig. 3). Therefore quantum entangle- 
ment is mainly manifest in the angular correlations. In 
fact, the nature of quantum entanglement can be appre- 
ciated from the probability distribution associated with 
the photon amplitude (Fig. 4). Along the k = q mani- 
fold, we see a narrow peak emerge near A9 = 9k—9q = TT. 
For smaller values ofba± (i.e., higher entanglement), the 
more pronounced is the peak observed. In other words 
kj^ is almost locked near —q± for sufficiently large trans- 
verse wavevectors. This demonstrates an interesting con- 
nection between quantum entanglement and localization 
in the transverse momentum space. In fact, recent stud- 
ies have begun to examine the role of entanglement in 
localization problems [2(1 l2l| . 

Fig. 4 suggests that wave vectors with large transverse 
magnitude are 'more entangled'. Therefore higher entan- 
glement may be achieved by simply choosing wavevectors 
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FIG. 4: A plot of the two photon amphtude square 
jCs(kx,q±)|^ in the k± = q± subspace. For convenience, 
the peak value is normalized to one. A cutting plane is shown 
for the entanglement enhancement scheme discussed in the 
text. The parameter is ba± — 0.25. 



with sufficiently large transverse components. This may 
be realized by adopting certain filtering procedures. As 
a demonstration, we consider an amplitude 

^(ki, qx) = a(kx, q±)0(fci - tic)0{q± - A*c) (7) 

in which low transverse wavevectors smaller than fj,c are 
removed by the unit step functions 6{x). We have per- 
formed a Schmidt decomposition of ^(k^,q^) numeri- 
cally. For the case ba± = 0.25 and /ic = 2a±, the corre- 
sponding Schmidt number is found to be if « 17.2. This 
is substantially higher than the original value K — 10.2. 
The enhancement is more drastic for the almost disen- 
tangled case ba± — 1. In this case we find iiT « 26 by 
using fj,c — cr±. Such an enhancement via amplitude fil- 
tering exploits the long tail of the sine function, and it 
has a cost in terms of detection rate. Physically it means 
that there is uncertainty of projecting the two photon 
amplitude onto ^(k_L,qj_). For the bcr± = 1 example 
above, the probability of finding the system in vA(k_L , q±) 
is about 6%. 

To summarize, we described a method to analyze and 
interpret physically realistic cases of high transverse en- 
tanglement in parametric down conversion, based on the 
Schmidt decomposition method. For the two classes of 
two-photon amplitudes defined by Cg and Cs, the de- 
gree of entanglement is controlled by a single parameter 
baj_ . We determined the Schmidt number as a function of 
ba± and identified two experimentally accessible regimes 
ba± 3> 1 and ba± ^ 1 in which values of entanglement 
higher than reported to date can be studied. In addition 
we proposed a method to enhance entanglement via large 
transverse wave vector selection. In fact one sees many 
possibilities to engineer entanglement through mode se- 
lection and manipulation of the shape of the pump pulse 
proffie The Schmidt decomposition presented here 
provides a basic framework for investigation of such pos- 



sibilities. 
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